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Abstract 

We have analyzed theoretically the photoemission spectra of the insulating com- 
pound Li2VOSi04. Recently, this compound has been proposed as the first exper- 
imental realization of the frustrated Ji — J2 Heisenberg model. Although it is well 
known that Li2VOSi04 is magnetically ordered in a collinear arrangement below 
T/v = 2.8K, there is some controversy about the coexistence of two collinear phases 
above T/v- Using a generalized t — J model we have obtained a complex spectral 
structure that can be traced back to the underlying collinear magnetic structures. 
We discuss the possibility to use ARPES experiments as a way to discern among 
the different scenarios proposed in the literature. 
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The Heisenberg model on a square lattice with antiferromagnetic (AF) ex- 
change interactions to first (Ji) and second nearest neighbors (J2) is the so- 
called J1 — J2 model. It has all the necessary ingredients to show a rich T = 
phase diagram: frustrating exchange interactions, low spin S=l/2, and low 
dimensionality. For J2/J1 > 0.55 the ground state is a twofold degenerate 
collinear phase (see Fig. [T]), with nearest neighboring spins aligned ferromag- 
netically along the x direction and antiferromagnetically along the y direction, 
or vice versa, characterized by the magnetic wave vector Q = (0, vr) or (7r,0). 
So, the 90° rotational symmetry of the square lattice is broken. Even if in 2D 
models with continous symmetry there cannot be a phase transition at finite 
temperature, Chandra, Coleman and Larkin (CCL)[2J argued that the twofold 
degeneracy of the collinear ground state could generate a finite temperature 
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Ising-like phase transition, since now there is a discrete symmetry that can 
be broken, although preserving the spin rotational symmetry. Recent series 
expansion calculations [3J, however, do not find evidence for such a finite T 
transition. On the other hand, it has been proposed [1] another mechanism 
in which the coupling of spins with the lattice is taken into account and a 
finite temperature Peierls-like transition may remove the twofold magnetic 
degeneracy. 
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Fig. 1. Collinear ordered ground states of the J1 — J2 model for J2/J1 > 0.55. 

These theoretical issues can be explored in the quasi-2D vanadium oxide 
Li2VOSi04, which has been found to be good experimental realizations of 
the J1 — J2 model[5]. In Li2VOSiOi, the S = 1/2 can be assumed local- 
ized at V^^{3d^) ions with competing superexchange interactions to first and 
second neighbors on a square lattice. NMR[5], neutron scattering [6j, reso- 
nant x-ray scattering |6], and magnetization measurements [5] indicate that a 
collinear structure is established below T/v = 2.8K. It should be noticed that 
the collinear transition at Tjy is triggered by the diverging 2D spin correla- 
tions present in the vanadium layers. Consequently, the J1 — J2 model is the 
appropriate one for the analysis of the magnetic properties of Li2VOSi04[5\, 
though the finite Neel temperature is due to 3D residual interactions. Above 
T/v the situation is less clear. It has been suggested by NMR experiments [5] 
that just above the Neel temperature the degeneracy between the two possible 
collinear phases is relieved by a structural distorsion. Furthermore, it has been 
pointed out[7J that magnetic domains of both kind of collinear phases coexist 
for T/v < T < Ji + J2. This coexistence has been deduced quite indirectly from 
the very low spin dynamics observed in NMR measurements. However, neu- 
tron scattering difraction experiments [H] do not reveal any structural phase 
transition. So, one question that remains to be answered is whether the rup- 
ture of the discrete degeneracy between the two collinear phases occurs at 
T/v or at a higher temperature T ~ Ji + J2, driven by a purely magnetic 
mechanism [2]. 

Consequently, it is worth to look for an alternative experiment able to detect 
different magnetic structures. For this reason, we propose that angle-resolved 
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photoemission experiments (ARPES), in which the induced photohole probes 
the underlying magnetic structure, can give relevant information that can be 
contrasted with the already existing experimental data. In general, photoemis- 
sion experiments are not used to extract magnetic information of a compound 
since the spectrum is not able to distinguish between short and true long range 
order. However, it has been recently observed [9] that the photohole propagates 
at different energy scales that are closely related to the underlying magnetic 
correlations, so that ARPES experiments could be used to discern among 
different magnetic structures, even short ranged ones. In particular, for mag- 
netic structures where ferro and AF correlations coexist, we have found [9] a 
many-body quasiparticle excitation that results from the coherent coupling of 
the hole with the AF magnon excitations, while the ferromagnetic component 
gives rise to a free-like hole motion at higher energy. 

In the present article, we compute hole spectral functions using a generalized 
t—J model appropriate for Li2VOSiOi. Using realistic parameters|8j we have 
found that the photoemission spectra clearly show different features when both 
collinear phases coexist, or not. This result could be helpful in elucidating 
the controversy about the magnetic properties of Li2VOSi04, in the regime 
Tn<T <Ji + J2. 



To compute the spectral function corresponding to a photohole induced in a 
photoemission experiment for Li2VOSi04,, we use a generalized t — J model 
that naturally takes into account the dynamics of a hole injected in the J1—J2 
model. This generalized t — J model can be written as 



H=-ti "^{cl^Cj^a + h.c.) - t2'^{cl^Ck,a + h.c.) + 
+Ji X! ^« ■ + +"^2 X! ^« ■ ^^^5 



where ti (^2) and Ji (J2) correspond to the hopping and the exchange inte- 
grals, respectively, to first (second) neighbors. The importance of hopping and 
exchange interactions to second neighbors is evidenced by the existence of a 
collinear magnetic order below and it is theoretically supported by LDA 
calculations [H]. In order to calculate the spectral function we use the spin- 
less fermion representation [10]. This is a standard procedure that leads to an 
effective Hamiltonian that explicitly includes the coupling between the hole 
and the magnon excitations of the frustrated magnetic order. The effective 
Hamiltonian can be written as: 
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where the operator hi. {a]j creates a hole (magnon) with momentum k. 
ek = 2ti cos ky is the bare hole energy, Mkk' = /^k'f k-k' — /^k^k-k' is the 
hole-magnon vertex interaction, where /3k = 2tisinA;a; + 4^2 sin fca, cos /cy, and 
Mk, Uk are the usual Bogoliubov coefficients, c^k is the spin wave dispersion 
in a coUinear order [2]. The expressions for Ck and /?k have been given for the 
collinear order (vr, 0), similar expressions hold for (0,7r). Treating Hef within 
the self-consistent Born approximation (SCBA)[10j, a self-consistent equation 
for the self energy Sk(ci^), at T = 0, can be derived 

Sk(a;) = ^ ^ I Mkq 1^ G'k-q(^ - UJq), 
q 

where the Green function is defined as Gk(t^) = (t^ — Ck — ^ki^))^^- Once the 
self energy is calculated, the spectral function can be computed as Ak(co') = 
-i/m GkH. 

Theoretically, at finite low temperatures, T < Ji + J2, it is expected a con- 
siderable 2D magnetic correlation length C, ~ exp(27rps/7'), where ps is the 
spin stiffness. Therefore, the low energy excitations at finite temperatures still 
resemble the magnon excitations of the ordered ground state. As the spec- 
tral function is insensitive to the long range order, a T = calculation that 
considers the coupling between the photohole and magnons will be valid for 
Tn <T <Ji + J2. 

Regarding the coexistence of both collinear phases in the range Tjv < T < 
Ji + J2, it has been suggested [7] that the size of the magnetic domains is pro- 
portional to the correlation length ^. If we assume that the domain sizes are 
sufficiently large, it is a reasonable approximation to consider that the photo- 
hole dynamics is only affected by the magnetic domain where it was induced. 
As a consequence, the expected spectral function in the case of coexistence of 
both orders will be the incoherent sum of the contribution of each order, that 
is, v4k(co') = -|- 74j°'''''(ci;)], as we assume that the photoelectron is 

extracted with equal probability from both kind of magnetic domains. 

To obtain the spectroscopic fingerprints corresponding to the coexistence, or 
not, of different magnetic domains for T/v < T < Ji + J2, we have evaluated 
the spectral functions for each case using the LDA estimated microscopic 
parameters for Li2VOSiOi: ti = —S.bmeV, t2 = —29.1meV, Ji = 0.83K 
and J2 = 9.81K^. In the upper pannel of Figure [2] it is shown the spectral 
function for a hole with momentum k = (vr, 0) propagating in the collinear 
order characterized by the wave vector Q = (tt, 0). This spectral function 
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has a finite lifetime resonance located near uj = 2ti, that can be associated 
with the hole propagation along the ferromagnetic chains in the y direction. 
At the bottom of the spectral function there is a weak quasiparticle peak, 
whose origin can be traced back to the coherent coupling of the hole with the 
magnons [T0f9] . We have observed that the broad resonance disperses along 
the ky direction hke e^, with a bandwidth ~ 4ti, while the quasiparticle peak 
exhibits a much smaller dispersion of order 2J2. 
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Fig. 2. Single-hole spectral function A(k, w) for k = (7r,0) and realistic parameters 
for Li2VOSiOi (see text). Upper panel: spectral function for the hole propagating 
in the (vr, 0) collinear order. Lower panel: spectral function for the hole propagating 
in a magnetic background of coexistent (vr, 0) and (0, vr) orders. 

As we mentioned above, for the case in which (vr, 0) and (0, vr) orders coex- 
ist (bottom panel of Fig. [2]) we have computed the spectral function as the 
incoherent sum of the contribution of each order. Due to the presence of the 
two kinds of collinear arrangements, the selected photohole with momentum 
k = (vr, 0), propagates preferently at two different energy scales, separated by 
an energy ~ 6ti, and which are related to the propagation along the ferromag- 
netic chains of each magnetic background. The differences between the spectral 
function for (7r,0), A(k, 0;)'^'^'°'', and the incoherent sum are most clearly ob- 
servable for momenta close to (vr,0) or (0,7r), while along the Brillouin zone 
diagonal, V (vr, vr), the spectral functions are equal for both orders. It should 
be noticed that, in the case of coexistence of phases, the discrete broken lat- 
tice symmetry is restored and, consequently, the spectral functions along the 
high-symmetry directions F (vr, 0) and F (0, vr) coincide. 

Another problem that still remains unsolved in Li2VOSi04 is the lack of an 
accurate estimation of the degree of frustration J2/ Ji- While there is a general 
consensus that Ji + J2 ~ 8.2i^, there is no agreement about the ratio J^jJx-, 
with proposed values ranging from 1.1 [5] to 12 [8]. Contrary to what happens 
with the magnetic susceptibility and the specific heat[8], we have observed that 
several observable features of the spectra are sensitive to the Hamiltonian pa- 
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rameters. For instance, within our scheme, the first moment of the spectrum 
for momentum k, that is, its center of gravity < e >k= / ix'A(k, uj)duj, is equal 
to the hole free hopping energy ek, with a bandwidth W = Ati. The second mo- 
ment of the spectrum < Ae >k= \J < e >kyA{k,u))du} characterizes its 
energy extension and it is equal to the average hole-magnon vertex interaction 
;^Z]q|-^k,qP ~ 3ti + 5^2- Finally, we have found that the weak quasiparti- 
cle signal at the threshold of the spectra disperses with a bandwidth ~ 2J2, 
because of the predominance of AF correlations between second nearest neigh- 
bors. As we can see, all the above spectral features are intimately related to 
the presence of collinear magnetic order, and, in principle, they can be read 
off directly from ARPES experiments, allowing an experimental estimation of 
the microscopic parameters. 

In conclusion, we have presented a theoretical calculation of the photoemission 
spectra for Li2VOSi04^, the first experimental realization of the J1—J2 model. 
We have used a generalized t—J model that takes proper account of the pho- 
tohole dynamics in a collinear magnetic background. We have shown that it is 
possible to extract spectral fingerprints of the underlying magnetic structures, 
directly associated with the ferro and AF components of the magnetic order. 
Therefore, we suggest that ARPES experiments would give further insight 
into open questions about this compound, complementing already available 
data. Of particular theoretical interest is the coexistence or not of the two 
collinear orders, (vr, 0) and (0, vr), in the temperature range T/v < T < Ji + J2. 
Our results point out that ARPES spectra can clearly distinguish between 
both situations, because the presence and dispersion of broad resonances in 
the spectral functions depend crucially on the underlying magnetic structures. 
We have shown that these broad resonances correspond to the photohole prop- 
agation along the ferromagnetic chains of the collinear domains. Moreover, we 
have found that several spectral features, like the quasiparticle and the cen- 
ter of gravity dispersions, would allow an experimental determination of the 
hopping and exchange energies. 
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